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Motivated by recent experiments on high-Tc cuprate superconductors pointing towards intra-unit- 
cell order in the pseudogap phase, we investigate three distinct intra-unit-cell-ordering possibilities: 
nematic, nematic-spin-nematic, and current-loop order. The first two are Fermi-surface instabilities 
involving a spontaneous charge and magnetization imbalance between the two oxygen sites in the 
unit cell, respectively, while the third describes circulating currents within the unit cell. We analyze 
the three-band Emery model of a single Cu02 layer including various on-site and nearest-neighbor 
interactions within a self-consistent mean-field approach. We show how these on-site and further- 
neighbor repulsions suppress or enhance particular lUC orders. In particular, we show that the 
attractive interactions necessary for nematic and nematic-spin-nematic orders in one-band models 
have their natural microscopic origin in the O-O on-site and nearest-neighbor repulsions in the 
three-band model. Finally, we find that while the nematic and nematic-spin-nematic orders cannot 
coexist in this framework, the loop-current order can coexist with nematic order. 

PACS numbers: 74.72.Kf, 73.22.Gk, 75.25.Dk 



I. INTRODUCTION 

Experimental evidence for various types of symme- 
try breaking in the pseudogap region of the phase dia- 
gram of the high-Tc cuprate superconductors has been 
accumulating in recent years. Neutron scattering experi- 
ments discovered a subtle staggered magnetic order in the 
pseudo-gap region of YBCOi and Hg-compounds^ that 
could be accounted for by either so-called nematic-spin- 
nematic orderi^ or circulating current loops^. On the 
other hand, neutron scattering^ and Nernst effect^ mea- 
surements on YBCO as well as SI-STM on BSCCO^ point 
towards an electronic nematic state. All these states re- 
tain the translational symmetry of the underlying crystal 
and can thus naturally be described by breaking intra- 
unit-cell (lUC) symmetries. Hence, identifying mecha- 
nisms for these symmetry-breaking possibilities and un- 
derstanding their competition is crucial for understand- 
ing the nature of the pseudogap phase. 

Theoretical investigations of translationally-invariant 
lUC order have so far been focussing on one partic- 
ular ordering at a time within simplified models each 
aimed at the ordering of interest. Nematic and nematic- 
spin-nematic order have only been studied in one- 
band models^i^i^rJ^ or in the extreme limit of infinite 
interactions^^. Loop currents, being more dependent on 
an lUC picture, have been studied in a mean-field picture 
with additional assumptions^ or numerically on small 
clusters or laddersi^r— . 

Here, we aim at a comprehensive investigation of 
lUC-ordering possibilities suggested by recent experi- 
ments.^ A theoretical description of intra-unit-cell or- 
der should be based on a microscopic model that al- 
lows for structures within the unit-cell. We thus start 
from a three-band model for the Cu02 plane, the so- 
called Emery model^, and consider various on-site and 
nearest-neighbor (nn) interactions (see Fig. [T]). We an- 
alyze three distinct lUC orders: nematic, nematic-spin- 



nematic, and loop currents. These phases can be distin- 
guished by the respective symmetries they break, both of 
the point group D^h and time reversal, as is summarized 
in Tab. ID For simplicity, only fourfold rotations, inver- 
sion, time reversal and combinations of these are shown. 
Within a self-consistent mean-field scheme we analyze 
and compare the origins of these phases and compare 
the infiuence of the different model parameters on them. 
In addition, we show how the 0-0 on-site and nn interac- 
tions result in effective interactions of (1^2 _y2 symmetry 
in one-band models, thus naturally leading to nematic or 
nematic-spin-nematic order. 

This paper is organized as follows: After introducing 
the Emery model describing the Cu02 plane, section Hill 
deals with nematic and nematic-spin-nematic lUC or- 
der through a decoupling of various interactions in the 
Hartree channel within self-consistent mean-field theory. 
In addition, we compare the three-band with the one- 
band model by focussing on the (partially filled) lowest 
of the three bands. Section IIVI examines lUC loop cur- 
rents by decoupling the nearest-neighbor interactions of 
the full Emery model in the Fock channel. Finally, sec- 
tion |V] compares the results for the different orderings 
and concludes. 
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TABLE I: The broken symmetries distinguishing the differ- 
ent lUC orderings with x denoting symmetries broken in the 
respective phase. For simplicity, we restrict the table to the 
fourfold rotation Ca, the inversion T, time-reversal operation 
T as well as combinations thereof. 
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II. MODEL 

The kinetic part of the Emery model^^ describing hop- 
ping of holes in the Cu02 plane is 

i.s y 

-^PpY. ^P\+0/2,sP^+0'/2,s+^■''■) 
i,s {u.v') 

i.s i^s u 

with and tpp the Cu-0 and 0-0 hopping integrals. 
Here, S}^ ^ creates a hole in the copper d^-^-y-i orbital at 
site i with spin s, v^^j^^j^ ^ creates a hole in the oxygen 

•py orbital at the site i + v/2 for i/ = x,y, and fif ^ 
are the corresponding number operators. The Cu sites 
i form a square lattice with unit vectors x and y, and 
the total number of lattice sites is N. The chemical po- 
tential ^ and the charge transfer energy A control the 
total and relative Cu/0 hole densities, and point 
to neighboring oxygen sites. 

In addition, we consider the interaction Hamiltonian 
including on-site interactions with strengths Ud and Up 
as well as nn interactions, Vpd and Vpp, 

i,u s,s' 

+^ppEEE"W.,"r+./2',.'- (2) 

The different orbitals and parameters of the model are 
shown in Fig. [TJ Setting tpd = 1, we fix the energy scale 
in the following. 




FIG. 1: The unit cell of the Cu02 plane with the copper 
d^2_y2 in the middle surrounded by the oxygen p^c and py 
orbitals. Also shown are the different hopping as well as in- 
teraction parameters used in the Emery model. 



III. NEMATIC AND NEMATIC-SPIN-NEMATIC 
ORDER 

For the above introduced Emery model, only the strong 
coupling limit, taking all interactions to infinity, has been 
analyzed for nematicity. Most theoretical investigations 
of nematic and nematic-spin-nematic order start from a 
single-band model, where in the weak-coupling limit a 
quadrupola r^'^'^" or a forward-scattering interactio n^ ^i^^ 
is introduced. For systems with a sufficiently high den- 
sity of states at the Fermi energy, e.g. due to a van Hove 
singularity, this can lead to a Pomeranchuck instability 
in the d-wave channel. Other studies of the Hubbard 
model without any additional (long-range) interactions 
found a nematic instability within a (weak-coupling) RG 
approacbi^, while DMFT calculations showed that the 
model maintains C4 symmetry^ii^ but becomes very 
susceptible to weak nematic driving fields (such as lat- 
tice distortions) close to the Mott transitions^ 

In the three-band model, the oxygen-oxygen nn inter- 
action prefers an imbalance in the hole densities of the 
neighboring oxygen sites, whereas the oxygen on-site in- 
teraction prefers to spin polarize the oxygen sites. The 
former interaction can thus lead to nematic order, the 
breaking of C4 symmetry, and the latter to either an 
overall magnetization on the oxygen sites or a nematic- 
spin-nematic order, which is invariant under a combi- 
nation of a C4 rotation and time reversal (see Tab. HI. 
Solving sclf-consistently the mean-field equations, we an- 
alyze the nematic and nematic-spin-nematic ordering in 
the following. 

A. Mean-Field Theory 

In this section, we focus only on symmetry breaking 
associated with the hole densities on the oxygen sites 
"^Zs ~ s)- -'^'^ absence of an overall magneti- 

zation on the oxygen sites, i.e. 

there are two distinct ways to break the lattice symmetry 
within each unit-cell, a nematic order and a nematic-spin 
nematic order. lUC nematic order in this model can be 
measured in terms of a spontaneous imbalance between 
the density of holes at the oxygen sites 

^^Kt + <i)-Kt + -^i)' (4) 

while nematic-spin-nematic order corresponds to having 
equal, but opposite magnetization on the two oxygen 
sites, 

Fig- m shows the distorted Fermi surfaces associated with 
these ordering possibilities as well as for ferromagnetic 
ordering for contrast. 
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FIG. 2: Static Fermi surface instabilities analyzed in this 
work: (a) Nematic phase breaking C4 symmetry, (b) nematic- 
spin-nematic and (c) ferromagnetic instability. In (b) and (c) , 
the solid and dashed lines denote the up- and down-spin band. 



Wc arrive in the following at the mean-field Hamil- 
tonian for nematic or nematic-spin-nematic order by a 
standard decoupling of all the interaction terms in Eq. ([2]) 
in the Hartree channel. We seek self-consistent solution 
with 77 ^ or 7^ 0, respectively, for nematic and 
nematic-spin-nematic order. 



1. Nematic order 

The mean- field Hamiltonian can be written in momen- 
tum space as 



(6) 



with = (plks'PLs'4j 



/ L 72(k) 7i(fc.)^ 

Hks = 72 (k) Jl{ky) 



and 



N 



-u, 



^PP~^ 



(7) 



(8) 



Here, n is the total density of holes, is the total density 
of holes on the oxygen sites, i.e., 

and the nematic order parameter 77 is defined in Eq.Q. 
The elements of the matrix ([7]) are given by 



72 (k) 



'2tpd cos—, 
-Atpp cos ^ cos 



and 



= A + Up— - Vpp- - /i, 

^y = A 4- ?7p— -t- Vpp- - /i, 

~ {n-nP) 
U = Ud 



(10) 

(11) 

(12) 
(13) 
(14) 



and also, we introduced the effective interaction param- 
eters 



Up = Up 



8Vpp 



SVr. 



pp — ^ y pp ^ pi 

Ud = Ud-Wpd. 



(15) 
(16) 
(17) 



In addition, we have put all the constant terms, i.e. 
2Vpdn — Vpd'n? , into the chemical potential fi. The mean- 
field Hamiltonian ^ can be diagonalized to yield three 
bands each with mixed orbital character and dispersion 
Caks, where a = 1,2,3 is the band index for the lowest 
lying and the two upper bands. 

In order to self-consistently determine the above intro- 
duced mean fields, we look at the grand potential per 
lattice site 



T 



log[l + exp(-eak.s/r)] 



N 



(18) 



For given values of -nP and ry, the chemical potential /x is 
implicitly given by solving 



dn N 

Q:,k,S 

distribution function n-p{x) 



(19) 



with the Fermi 
l/(exp(x/T) -I- 1). Self-consistency equations for and 
1] are found by extremizing the grand potential to 
be 



N{Up + 2Ud) : 



k,; 



"f(Co 



2Udn 



dnP 



U„ 



and 



'lF(Caks)- 



PP 



a,k,^ 



dr] 



2Ud 
(20) 

(21) 



Second-order phase boundaries for nematic ordering 
can be determined using the stability condition of uj by 
requiring 



drf 



( Ypp _ i_ 
V 4 N ^ 



J7=0 



1 



\ dn ) 



. 4Tcosh^%^ 

k.s 2T 



drf 



dr] 



0. (22) 



i)=0 



This is equivalent to analyzing the linearized self- 
consistency equation for rj. To additionally find first- 
order phase boundaries requires examining the grand po- 
tential for the global minimum. 



2. Nematic-spin-nematic order 

In complete analogy to the mean-field decoupling intro- 
duced above for the nematic order, we find the mean-field 
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FIG. 3: The critical strengtli of the oxygen-oxygen interaction 
Vpp needed in order to enter a nematic phase as a function 
of hole density n for tpp = 0. For numerical reasons, the 
calculation has been carried out at T = 5 • 10^^[tpd]. The 
dashed lines denote the values of Vpp used for Fig. [S] 



Hamiltonian for the nematic-spin-ncmatic ordering with 
non- vanishing r]s as defined in Eq. ([5]). We only have to 
replace the diagonal elements of the Hamiltonian ([7]) by 
the (now spin-dependent) 



(23) 
(24) 
(25) 



Up 


T 


-sUp- 


Up 


T 


+ sUp^ 
^ 4 



-- {n - nP) 

(ids — Ud /i, 



and 



Up^-Ud- 



(26) 



The interaction parameters Up and Ud arc again given by 
Eqs. p5|) and ([T7|. respectively. Note that now, the in- 
teraction driving the instability is not the oxygen-oxygen 
nn interaction Vpp, but the oxygen on-site interaction Up. 
This ordering is thus in direct competition with an overall 
magnetization m on the oxygen sites as given in Eq. Q . 
The critical Up for a ncmatic-spin-nematic instability to 
occur is again determined by analyzing the stability con- 
dition for the corresponding grand potential, d^^to = 0. 



B. Results 

1. Nematic order 

Our goal is to investigate the effect each parameter has 
on the nematic instability. For this, we use the linearized 
self-consistency equation p2|) to map out various phase 
boundaries in the parameter space. Due to the large 
parameter space of the three-band model we present re- 



FIG. 4: Critical interaction strength for different values of the 
oxygen-oxygen hopping tpp. 



values for the cuprates. Realistic values for the 0-0 hop- 
ping and the Cu on-site interaction are tpp « 0.2-0.5 and 
Ud ~ 6-8i2i*2^ After calculating a general phase diagram 
and looking at the influence of finite 0-0 hopping, we 
can thus for example analyze the effect of the Cu on-site 
interaction on the nematic phase formation. 

lUC nematic ordering within this mean-field theory 
arises through a Stoner-type instability. It therefore re- 
quires a finite interaction strength for all hole densities 
away from the van Hove filling, where the diverging den- 
sity of states allows for an instability at infinitesimal 
Vpp. Fig. [3] shows the critical oxygen-oxygen interaction 
strength V^^ needed to enter a nematic phase at T « 
for Cu on-site interaction strength Ud = ^ and tpp = 
with the van Hove singularity at rivH = 1- Note that the 
doping dependence of V^^ is not symmetric around the 
van Hove point. This is a multi-band effect and we will 
return to this in section UlI CI 

The main effect of a finite oxygen-oxygen hopping tpp is 
to shift the van Hove singularity to higher hole densities 
as can be seen in Fig. 31 where we show the hole-density 
dependence of V^^ for different tpp. In addition we see 



suits with A = 2.5, U„ 



3 and Vm 



1 as realistic 



that the nematicity is gradually suppressed upon an in- 
crease in tpp, which reflects the fact that a finite tpp adds 
to the two-dimensionality of the system. For the rest of 
this section, we will focus on the case of tpp = 0^ 

We now turn to the T-n phase diagram shown in 
Fig. [5ja). For the phase diagram, we obtain the second- 
order phase boundary from the linearized self-consistency 
equation for Vpp = 2, 1.75 and 1.5 (dashed horizontal 
lines in Fig. For Vpp ~ 1.5, a small dome almost 
symmetric around the van Hove filling is found while for 
higher Vpp, the dome becomes asymmetric with respect 
to the van Hove filling rivu = 1- At low temperature, 
we expect the phase transition to be of first order as was 
shown for one-band models within mean-field theory in 
Refii^. Explicitly examining the full grand potential for 
Vpp, we indeed find second-order transitions at higher 
temperature, i.e. for hole densities in the middle of the 
dome and first-order transitions for densities at the bor- 
der of the dome, as indicated for Vpp = 2 by the solid 
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FIG. 6: Phase diagram for Vpp = 1.75 and difTerent Cu on- 
site interaction strengths Ud ~ 8, 9, 10. Shown are again only 
'second-order' phase boundaries. 



of Figs. [Sfa) and [HI As discussed in section IIII CI this 
is due to an increased hole density at the oxygen sites 
as well as a reduced level separation between the lowest 
lying bands for larger Ud- 



2. Nematic-spin-nematic order 



FIG. 5: (a) Phase diagram for the different values of the O- 
O nearest-neighbor interaction Vpp = 2, 1.75 and 1.5. At 
low temperature, there would be first-order transitions, only 
shown for Vpp = 2 by the solid lines, before the normal state 
becomes unstable (dashed lines). Figs, (b) - (d) show the 
free energy as a function of 77 for n = 1.05, n — 1.095 and 
n = 1.105 at T = 0.001, illustrating the first-order character 
of the low temperature transition, (e) The nematic order 
parameter 77 as a function of temperature for Ud = 9 and 
n = 1.05 showing the second-order transition at Tc = 0.027. 

lines in Fig. [5ja). Note that in the case of a one-band 
description it was shown that fluctuations may make the 
first-order transitions continuous^^ and we thus expect a 
similar behaviour here. 

To illustrate the first-order character of the low- 
temperature transition, Figs. [5jb)-(d) show the free en- 
ergy as a function of ?; at T « for a hole density n deep 
inside the phase, where the normal state is metastable, 
and where the nematic state is metastable. The second- 
order character of the transition on the top of the dome 
is best seen in the T dependence of the order parameter 
77. This dependence is shown in Fig. [5jc). 

What is particularly noteworthy from our survey of pa- 
rameter space is that the Cu on-site interaction Ud tends 
to stabilizes the nematic phase as shown in Fig. [51 In 
fact we find that the effect of increasing Ud is almost the 
same as increasing Vpp as is apparent upon comparison 



We only present the doping dependence of the critical 
oxygen on-site interactions Up which drives the nematic- 
spin-nematic (see Fig. [7]) , as technical details of the self- 
consistency analysis for the nematic-spin-nematic phase 
are very much analogous to that for the nematic phase. 
For completeness, we also show the critical Up for a fer- 
romagnetic instability with m as defined in Eq. ([3]) and 
i]s = 0. This competing ferromagnetic instability is only 
favored over the nematic-spin-nematic order for n < 1 
(electron doping). Notice that the magnitude of Up for 
the nematic-spin-nematic order to occur is almost an or- 
der of magnitude larger than the magnitude of Vpp for 
the nematic instability. This is due to the fact that here, 
the holes on an oxygen site with spin s are only interact- 
ing with the holes on the same site with opposite spin, 
while for the nematic instability, the holes interact with 
the total hole density of all four neighboring oxygen sites. 

Now we can compare the influence of various interac- 
tion strengths for nematic-spin-nematic ordering to that 
for nematic ordering. We find that Ud, A and tpp have 
the same effect for both types of ordering: Ud increas- 
ing and A and tpp decreasing the tendency towards both 
orders. However, increase in Vpp leads to larger U^ as it 
reduces hole occupation of the oxygen sites. We thus find 
that the interaction driving the nematic or nematic-spin- 
nematic instability hurts the occurrence of the respective 
other phase. 
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FIG. 7: Critical oxygen interaction strength for a nematic- 
spin-nematic (rjs) and a magnetic (m) instability on the oxy- 
gen sites. Here, Ud = 9, Vpd = Vpp = 1 and tpp — 0. 



while the oxygen nn interaction only has a spin- 
symmetric part. We now start by expressing the spin- 
symmetric part of the inter-oxygen interactions in the 
basis Caks- The terms that only involve lowest bands are 



7L7L')^'k^'k''^iksnik', 



k.k' s.s' 



k.k' s.s' 



(32) 



with fiiks = ^iks^iks ^^^^ density operator for the lowest 



lying band. Finally Eq.([32]) can be reorganized into an 
isotropic part 



i^im^k^'k'^ik.nik'. 

k.k' s.s' 



(33) 



C. Comparison to one-band model 

In this subsection, we highlight similarities and differ- 
ences between the mean-field theory of the three-band 
Emery model and previous studies of effective one-band 
models. We first derive effective attractive interactions in 
spin-symmetric and antisymmetric channels with (i^2_j,2 
symmetry for the lowest band of non-interacting model, 
in terms of repulsive oxygen-oxygen interactions. We 
then discuss the multi-band effect in the "Stoner-like" 
self-consistency condition. 

In order to see how Vpp and Up lead to effective at- 
tractive interactions for the lowest lying band £,1^^, we 
express the interaction term in the diagonal basis Caks 
of the non-interacting Hamiltonian (setting tpp = 0). In 
this basis the oxygen operators read (see App. |A|) 



Pxks 
Pyks 

where Uk = 



-7l2:^^kClks 
-7la"kClks 



7ll/C2k,s + 7lx^kC3ks, 
7l2:C2ks + 7ly'"kC3ks, 



Wk 



sin ^ with 



Wk = arctan 



2V7i'(fc.)+7i'(fc,) 



(27) 
(28) 



(29) 



where Lip = {Up + 8Vpp)/2, and the effective "F2,s" part 



Wk^k'^lks^llk'. 



where V, 



pp 



k.k' s.s' 

(SVpp -Up)/ 2 and 

(cos kx — cos ky) 



dk = 



(2 -I- cos kx + cos ky) ' 



(34) 



(35) 



Some remarks are in order. Ea. p4p explicitly shows 
that repulsive Vpp leads to effective attractive interac- 
tion that can drive nematicity for the lowest lying band. 
It also shows that Up hinders nematic ordering. Further- 
more, we see that these inter-oxygen interactions are act- 
ing only on the portion of hole density in the lowest lying 
band ^^^.^ that can be attributed to oxygen occupation 
since the oxygen occupation number 



' = i^E-knF(aj- 

k,s 



(36) 



Following the same procedure, the spin-antisymmetric 
part of Eq. (|3ip can be organized into an isotropic part 
and the effective "i^2,a" part: 



and 



Hi 



iiih) 



Vliikx)+ji{ky)' 



x,y 



(30) 



with 7i(fci) = —2tpdCos^ as it was defined in Ea.([TU)>. 
Note that the oxygen on-site interaction can be separated 
into spin-symmetric and antisymmetric parts as 



kt'^i^k'i 



k.k' v=x,y 

u 



2N 



^ '^i'Tiiyksni^k's' - ss'niyksni,k's'), (31) 



k,k' i'=x,y s,s' 



^k^k'^lks^lk', 



k,k' s,s' 

U 



^^^ss'dkdu'vlvl,niusniu's'- (37) 



8N 

k.k' s.s' 

We therefore find explicitly that oxygen nearest-neighbor 
and on-site interactions in the three-band model lead in 
a one-band model to attractive spin-symmetric and an- 
tisymmetric interactions of dj.2_j^2 symmetry, i^2,s and 
F2.a, driving nematic and nematic-spin- nematic order, 
respectively. 

We now turn to the multi-band effect in the linearized 
self-consistency equation ([2^ . For the parameter space 



7 



a) 



b) 



G°'°'(k,iwn) G««(k,ia>„) 

{5H)c,c,'(^^^^^^^^^ {5H)aa {SH)af3 

FIG. 8: The two bubble diagrams involved in the linearized 
self-consistency equation, where we have used the short nota- 
tion {sn)c,fi = {a\sn\i3). 




of interest to the cuprates, only the lowest band of the 
(mean-field) Hamiltonian with energy ^iks is filled at low 
temperatures. Hence, Eq. ([2^ amounts to 



FIG. 9: The different current patterns arising from the opera- 
tors ^1-43 in Eqs. (gll) and (gSJ and ©1-43 in Eqs. (gg]) - (|48)) . 
Combining A2s with leads to the loop-current phase Qi, 
while A33 (.445) combined with D^s (©48) leads to On. 



( Yee. _ J_ 

V 4 N ^ 



4rcosh^ %f 



Iks 



lk.< 



19772 



0. (38) 



While the first line is the familiar result from one-band 
mean-field calculations with the second term being the 
familiar polarization bubble, the term on the second line 
has no analogue in simple single-band models. This term 
grows with total hole density and is thus responsible 
for the asymmetry around the van Hove filling found in 
Secinill 

To better understand Eq. ([38]), we interpret ne- 
matic order as a perturbation around the isotropic 
[f] ~ 0) Hamiltonian H^^f^ (as done in App. |X] for 
tpp = 0) and write the derivatives in Eq. (|38p in terms 
of the perturbation-theory expansion parameters, see 
Eqs. (|A13|) - (jAlSp . Introducing the (isotropic) Green's 
functions Gi"'°^(k,w„) = (iw„ - Cik°V^ "^'^^^ ^aks = 
S.a'ks\ri=o the eigcneuergies of ^^8°^ ^^'^ using the rela- 
tion 



^GL-°)(k,c.„)GL-°)(k-Kq, 
_ '^F(Cks )-"F(eik+qs) 



^(iso) _ ^(iso) 



q^O 



(39) 



we find 



A^f^ 4rcosh2%f V di] 



77=0 



= j^T.Gt\^,^n)Gf^°\\^,uj,mSn\l), (40) 

ks 

the familiar polarization bubble as depicted in Fig. |Sl[a) . 
For the case of multiple bands, also inter-band inter- 
actions should be taken into account as indicated by 



Fig. IHb). The corresponding expression yields 
J2 E Gl^'°HKu^n)Gf°\k,u;^)\{a\5nm\' 

\{a\SHm' 



a^/i k 



f (iso) _ f (iso) 
^aks ^/3ks 



1 x^^Ciks 



k.s 



, (41) 



where in the last step we have used Eq. (|A15p and the 
fact, that only the lowest band is occupied for low tem- 
peratures. The asymmetry is thus a multi-band effect 
unlike the asymmetry found in RefJ^, which is due to 
an asymmetric density of states. It is now clear why 
the asymmetry only appears for large values of Vpp: only 
when the 0-0 nn interaction is comparable to the band 
separation, the influence of this term becomes visible. 
Increasing Ud has then two effects, both enhancing ne- 
maticity: in addition to increasing the hole density on the 
oxygen sites, it shifts the lowest band slightly up in en- 
ergy, decreasing the level separation to the second band, 
hence increasing the importance of the inter-band term. 
This explains why an increase in Ud has such a similar 
effect as an increase in K.^- 



IV. LOOP CURRENTS 

In the loop-current picture, the staggered magnetic 
moments observed in experiments^ originate in circulat- 
ing electron currents around 0-Cu-O triangles within the 
unit cell. This phase was introduced and analyzed in a 
mean-field approach by Varma<^ Stipulating a Cu-0 hop- 
ping integral that depends on the hole doping and a van- 
ishing charge-transfer gap, a phase-diagram was drawn in 
qualitative agreement with the pseudogap phase. While 
later exact diagonalization calculation on small clusters 
concluded that the energy scale of such current loops was 
too small to account for the phenomena associated with 
the pseudogap rcgioufiS, RG calculations for two-leg lad- 



8 



ders found incommensurate loop currentsii and a vari- 
ational Monte Carlo study found that the 0// current 
pattern is stabilized in intermediate system sizea^^. 

In this section, our aim is not to answer the question 
whether such loop currents exist in the parameter range 
usually assumed for the cuprates, but again to investigate 
the microscopic origin and the influence of the various 
model parameters. This then allows for a comparison 
with the two ordcrings of the previous section. 



A. Mean-field theory of loop currents 

To analyze this loop phase, we perform a similar cal- 
culation as in Refi^, however with some important dif- 
ferences: First, we only use the Hamiltonian as defined 
in Eqs. ^ and ([5]) without any implicit assumption re- 
garding doping dependence of parameters. Further, our 
calculation includes the 0-0 nn interaction, which we 
decouple analogously to the Cu-0 nn interaction. Fi- 
nally, wc do not reformulate the Hamiltonian in terms of 
flux through the Cu-0 triangles, but decouple the nn in- 
teraction terms and directly solve for the self-consistent 
mean-field solution. Our conventional treatment of the 
problem leads to different doping dependence and an ad- 
ditional Vpp contribution compared to Ref. — . 

We start with the interactions in Eq. ([2]) and follow 
the Cu-O-interaction decoupling of Varmai^ by defining 
the operators 



A 



l/2is 



± idlsP,+y/2,s + dlsP^-y/2,s)\ ' (42) 



•^3/4i, 



(43) 



2 ('^L'5-Pi+x/2,s '^L'5Pi-£/2,s) 

± idi,sPi+y/2,S " di,sPi-y/2,S 

Introducing the (complex) mean-field order parameters 
R.e''^" = Vp^Y^iA^s) (44) 

s 

with (Aus) — (Aiyis) independent of site i, the Cu-0 in- 
teraction can straight-forwardly be mean-field decoupled 
in the Fock channel 

— Vpd A^igA^j^^, ~ 

iM s.s' 

~^(i?,e-*^A..+h.c.) + ^^i?2. (45) 

The order parameters i?^ correspond to the current pat- 
terns shown in Fig. |9] and can again be classified accord- 
ing to the symmetries they break. We first note that 
can not lead to a stationary current loop. Focussing for 



the other order parameters again only on C4, I, 7" and 
combinations thereof, we find that the order parameter 
arising from A2S, corresponding to 0/ in Refi^, differs 
from the order parameters arising from Azs and ,44s cor- 
responding to the 0/7 phase, in that it preserves C4 o 7" 
(see Tab. p. 

Next, we similarly look at the 0-0 interaction. For 
the decoupling, we again introduce operators of distinct 
symmetry, 



'^'l/2is ^{p'i-x/2,sPi+y/2,s ~^ Pi+x /2,sPi+y /2.s 



' Pi+x/2,sPi-y/2,s "F Pi-x /2,sP i-y /2,s 



and 



+x/2.sPi+y/2,s 



■Pi 



i-x/2,sPi-y/2,s J ' 



(46) 



(47) 



^Ls ~ '''{Pi~x/2,sPi+y/2.s P\+x/2,sPi-y/2,s]- (^8) 



This allows us to introduce (site-independent) mean- 
fields 



(49) 



and decouple the 0-0 interaction term in the Fock chan- 
nel as 



i^U S,S' 

1 ^ 



2 ^ 



h-C-) + ^E(^')'- (50) 



2K 



Looking at Fig. |9l we see that Vis only breaks time- 
reversal symmetry, while 7)28 leads to the 0/ phase and 
T^'i/is to the 0/7 phase with the respective broken sym- 
metries. 

In the following, we are only interested in the phase 
077 and therefore only keep the two order parameters 
R3 = R and R^ = Rp, which have the same symmetry 
and mix, finite. Analogously, we could also choose R4 and 
i?4 (see Fig. [9]). An order parameter yielding a current 
(rather than a bond-density) has to have an imaginary 
part and for simplicity, we set the phases to 03 = 03 = 
/2,2a Note that non-zero R and Rp, while corresponding 
to non-zero loop-currents, do not guarantee the absence 
of macroscopic currents. 

Contrary to Sec. IIIH where the mean-field decoupling 
led to a shift of the diagonal elements in the Hamiltonian, 
here it leads to new hopping elements in the Hamilto- 
nian d?!) 



71 (fci ) = - 2tp<i cos y - ii? sin ^ 



2 ' 



(51) 



and 
72 (k) 



-4:t„„ COS — cos — 
^^2 2 

■ 7-» / • • \ /rr»\ 

— i/t„(sm — COS — — cos — sm — I. (52) 
^^ 2 2 2 2 
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FIG. 10: Critical interactions (Vpp, V^i) for Ud = 9, f/p = 3, 
tpp = 0.1, A = 2.5 and different hole densities. The dashed 
and dotted lines for n = 0.9 illustrate the influence of the Cu 
on-site interaction and the charge transfer gap. 



Ud results in higher critical interaction strengths. How- 
ever, this mainly affects V^^ due to the change in the 
oxygen hole occupancy, while V^^: is almost unchanged. 

Contrary to the nematic and the ncmatic-spin-ncmatic 
order of the previous section, the critical interaction 
strengths here are monotonically decreasing with increas- 
ing hole density (see Fig. ITU)) . This is due to the fact that 
the current loop in a mean-field approach arises due to a 
Fock-type rather than Hartree-type decoupling and hence 
not a Stoner-type instability. As the whole dispersion is 
altered by the decoupling, increasing the hole density in 
the lowest band increases the tendency towards loop cur- 
rents. In order to find a phase diagram as found in the 
cuprates, additional assumptions to the model have to be 
made, such as a density-dependent hopping, e.g. of the 
form tpd — >■ 'tpd\x\ with x = n— 1, as in Varma's analysisi^ 



DISCUSSION AND CONCLUSIONS 



For simplicity, we decouple the nearest-neighbor interac- 
tion only in the loop-current (Fock) channel and thus the 
diagonal elements yield 

C,,y = A + C/p— -/i, (53) 

= £/ / ^ ^ -/i, (54) 

and 

f{nP,R,Rp) _ {nP f {n - n^f 

N 8 ^ Vpd 2Vpp' 

(55) 

To investigate the occurrence of instabilities, we need to 
account for the fact that the two order parameters R and 
Rp are coupled by symmetry. This means that we can 
not investigate their respective instabilities separately as 
done for the order parameters in the previous section. In- 
stead, the pairs of critical interaction strengths (V^^, V^^) 
are given by a vanishing eigenvalue of the Hessian matrix 
of the grand potential a;. 



B. Results 

The solid lines in Fig. [TOl show the curves (V^p, V^^) for 
Ud = 9, Up = 3, tpp = 0.1 and different hole densities. 
Due to coupling of the two order parameters R and Rp, 
the critical Cu-0 nn interaction V^^ is reduced by a finite 
Vpp. We again study the influence of different parame- 
ters on the critical interaction values. As the dashed and 
the dotted lines for n = 0.9 show, increasing the charge 
transfer gap A or reducing the copper on-site interaction 



Starting from a three-band model and applying a 
mean-field approach - despite its obvious shortcomings - 
we gained valuable insights about the microscopic repul- 
sive interactions that can promote various lUC orders. 
We found that the Cu on-site interaction Ud increases 
the tendency towards all the studied orderings by shifting 
more holes to the oxygens. The charge transfer gap A has 
the opposite effect. Also, different interaction parameters 
affect the different instabilities differently: while the O 
on-site repulsion Up only favors the nematic- spin- nematic 
phase and the Cu-0 repulsion Vpd the loop currents, the 
nearest-neighbor 0-0 repulsion Vpp helps both, the ne- 
matic and the loop-current phase (see Tab. [Ill for a sum- 
mary of all the model parameters). Further, we could 
microscopically motivate attractive interactions _F2,s and 
^2,a with a repulsive (longer-ranged) 0-0 repulsion Vpp 
and O on-site repulsion Up, respectively 

A comment on the magnitude of the interactions neces- 
sary found here is in order: the energy scale of the pseu- 
dogap phase in the cuprates is of order 100 Kelvin. For 
any of the above phases to reach to such high tempera- 
tures, unrealistically large interactions are needed within 
our mean-field calculation. Looking at the nematic phase 
for example and taking tpd ~ leV, Vpp needs to be of 
order 2tpd as can be deduced from Fig. [5]Ja). Also, in- 
creasing the Cu on-site interaction strength Ud to enter 
a nematic phase leads to unphysically large values. For 





Ud 


Up 


Vpd 


Vpp 


tpp 


A 


nematic 


+ 






+ 






nematic-spin-nematic 


+ 


+ 










0/1 loop current 


+ 




+ 


+ 







TABLE II: Summary of the effect of the different parameters 
in the Emery model on the different lUC orders, where -|- 
denotes a parameter that helps a specific order and a - denotes 
a hindering parameter. 
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the nematic-spin-nematic phase, the respective interac- 
tions need to be even larger as can be seen in Fig. [T] 
Finally, for the loop-current phase, we deduce values for 
the critical interaction strengths from Fig. [TUl that are 
much higher than realistically expected. However, the 
aim of our analysis is not to answer whether this phases 
exist in the cupratcs - a mean-field analysis would cer- 
tainly not be the appropriate approach for such a task 
- but to analyze this different lUC ordcrings within the 
same framework. 



We can also draw some conclusions about the compe- 
tition or coexistence of the lUC-ordered phases from our 
calculation. Even though symmetry allows the two orders 
discussed in Sec. IIIIl the nematic and the nematic-spin- 
nematic, to coexist independently, they are promoted 
by different interactions, Vpp and Up, each hurting the 
respective other phase. In addition, both depend in a 
mean-field picture on the presence of a van Hove singu- 
larity. Even if the interactions were tuned in a way as to 
allow for both instabilities, having one kind of order re- 
moves already the high density of states from the Fermi 
level, thus preventing the system from entering the other 
phase. In contrast, the loop-current phase is promoted by 
the same interaction as the nematic phase, Vpp, and docs 
not depend on a high density of states at the Fermi level. 
A deformation of the Fermi surface has thus no direct in- 
fluence on this instability. Being of different symmetry, 
a nematic and a loop-current phase can therefore coexist 
in general independently. 

The mean-field analysis and our exploration of the rich 
phase space of the three-band Emery model in this pa- 
per can serve as a stepping stone towards more sophis- 
ticated calculations of lUC orders and their interdepen- 
dence. For instance, extension of the calculations inl^J^ 
to the case of three bands might provide further valuable 
insight. In particular, investigation of the interplay be- 
tween these lUC orders and superconductivity in a gen- 
uinely strong coupling approach will be of great interest. 
As superconductivity cannot be accessed within mean- 
field theory with purely repulsive interactions, we left 
out this important issue altogether—. 

The possible coexistence of nematic and loop-current 
phases we find in this work is interesting in light of exper- 
imental observations of both lUC nematic orderi^ and 
lUC staggered magnetism in underdoped YBGO^. On 
the other hand, in Hg-compounds only lUC staggered 
magnetism has been observed^. In order to test whether 
coexistence of both orders is a generic feature, we pro- 
pose a measurement of anisotropy in Nernst effect in the 
presence of a symmetry breaking field on Hg-compounds. 
For example, an in-plane magnetic field could align pos- 
sibly existing nematic domains in Hg-compounds which 
is tetragonal otherwise, much in the same manner as in 
SraRuaOT^. 



Appendix A: tpp — analysis 

In this appendix, the case tpp = is investigated for 
?/ ^ for which we write the Hamiltonian ([7]) as 



(iso) 



ks5 



where 



n 



(iso) 



/ 7i(fc.) 

liiky) 



(Al) 



(A2) 



and 








(A3) 



with 



Sp 



A + Up— - fi. 



(A4) 
(A5) 



The unperturbed Hamiltonian (jA2p can straight- 
forwardly be diagonalized yielding the eigenenergies 
£^ = ' and 



1 



(iso) 

3/lks 2 ^"^P 

± 



\lli6''°^-^T°^r+lKk.)+lKky) (A6) 



with the corresponding states \va) given by the eigenvec- 
tors 



Vl 











( llxU-k \ 






Co I 




7ly'"k 










I «k / 



(A7) 



In these equations, we introduced Uk = cos^, Uk 



with 



Wk = arctan 



V7i'(fcx) + 7f(fcy) 



.^(iso) 



and 



llxiv) 



7l(fcx(y)) 



f(iso) 



^lt{kx)+lt{ky) 



(A8) 



(A9) 



For ry — 0, we can thus express the eigenenergies in pow- 
ers of ?7 in a text-book perturbation-theory expansion, 

Cok. = e + ni^2s + ^'dll + 0(ry3), (AlO) 
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with 



and 



^(2) 



\{vc.\sn\vp)\^ 



E 



^(iso) 
P^a Saks 



(All) 



(A12) 



The derivatives appearing in the self-consistency equa- 
tions ((20)) . ((2T|) and ([22]) can thus ah be expressed ana- 
lytically through Eq. (|A10p yielding 



(A13) 
(A14) 
(A15) 







_ '^?Qks 


dnP 




dnP 


dr] 


?)=0 




drf 


ri=0 


- 2f^'^ 



Evaluating the derivatives for the oxygen hole density in 
Eq. (HOI), we find using (|JT3| 



Iks 



dnP 



dnP 



dnP 



i)=0 



i)=0 



i)=0 



4 2 2 ' 



,C/p , Ud, Ud 



(A16) 
(A17) 
(A18) 



Using these derivatives, the self-consistency equation for 
the oxygen occupation number simplifies to 



k,s 



ISO } 

ks 



(dk:')} (A19) 
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